Lattice Models for Magnetic Fluids: Correlations Between Order Parameters 



Dorina ANDRU VANGHELlj], Horatiu COVLESCU 
Gheorghe ARDELEAN, Cristian STELIA 
Department of Theor. and. Comp. Physics, 
West University of Timi§oara, V. Pdrvan Str., no. 4, 
1900 Timisoara, Romania 

Abstract 

Magnetic fluids are colloidal suspensions of ferromagnetic particles covered with a surfactant 
layer, dispersed in a host liquid. The existence of cooperative phenomena in such magnetic colloidal 
systems, makes the determining of their physical properties a difficult task. In the framework of 
statistical mechanics, taking into consideration interparticle interactions, the study of equilibrium 
properties and phase transitions of magnetic fluids can be best done. 

In this paper we analyze some basic aspects about the real possibility of describing the critical 
properties of these systems in the language of spin-1 lattice-gas model and we also propose the use 
of a higher spin (3/2) model. 

The main features of magnetic fluids are related to their strong magnetic properties (magnetic 
susceptibility, \ ~ 0.1 - 1). 

For a long period of time, all theoretical approaches in the field relayed on hidro dynamical models, 
thus regarding magnetic fluids as continuous media. The absence of interparticle interactions and 
the postulation of the Langevin law as a magnetic state equation, common to all this models, let 
unexplained a number of phenomena. Within these, of great practical and theoretical importance is 
the aggregation process, i.e. the formation of large agglomerates of magnetic particles in the presence 
of a magnetic field. In most cases, these agglomerates vanish when the field is switched off, which 
suggests that the process is reversible. Meanwhile, deviations of magnetic properties (magnetization, 
initial susceptibility) from the Langevin behavior were observed as a consequence of aggregation. The 
phenomena was recorded by various experimental methods (see for example Q). 

A new approach, from the viewpoint of statistical mechanics, which takes into consideration the in- 
terparticle interactions, was needed. Thereafter, statistical models appeared in the theoretical studies 
on magnetic fluids. 

Cebers, |^], and, independently, Sano and Doi, ||, emphasized that the occurrence of agglomerates 
in the presence of an external magnetic field could be regarded as a phase separation generated by the 
thermodynamic instabilities of the colloidal solution. The description of these phenomena was given 
in the framework of the mean field theory. Other approaches include the cell model Q and the hard 
dipolare spheres model Q. These models led to remarkable results on the magnetic fluids equilibrium 
thermodynamic and on their phase transitions. 

Another general approach was made by using the lattice-Boltzmann model on a triangular bidi- 
mensional lattice, ||. In fact, te model is a particular case of the Ising-^ model. Is is well known, that 
for the Ising-^ triangular bidimensional structure, critical properties are determined by the fixed point 
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Hamiltonian, it's dipolar coupling constant being J = 0.34 in the absence of an external field. The 
density of particles for this value of the coupling constant is n = 10 24 ferroparticles/m 3 |7|, a result 
which is in good agreement with the critical value of the volume fraction. As a consequence of this 
observation, we can conclude that the Ising- 1 model can be successfully used to analyze, qualitatively, 
both critical behavior and stability properties of magnetic fluids. 

The development of theories on exact or approximate solutions of the Ising- 1 ; models is based on 
the existence of a single order parameter, identified with the density, concentration, crystalografic 
order parameter, magnetization or polarization of the system. 

But, in the case of magnetic fluids, a good description of the agglomeration phenomena can not 
be completely specified by a single order parameter, because, in this case we have to deal with two 
types of ordering processes: magnetic and structural. In these conditions, we will need at least two 
order parameters, which can be done in using a Ising-1 lattice model (||, 

If we consider the following general Ising-like Hamiltonian 

H = -Hj^Si - - D^(Si) 2 - Kj2(Si)\S,) 2 (1) 

i i,j i i,j 

where Qi = (Si) 2 is the singleparticle quadrupolar operator, as the two order parameters we shall use 
the quantities 

M =< Si > (2) 

and 

Q =< (S l ) 2 > (3) 

These two parameters are not independent in the sense that, generally speaking, M ^ implies Q ^ 0. 

In the Bragg- Williams mean-field approximation, for the spin-1 general Hamiltonian, where H = 0, 
the kinematical coupling between M and Q is determined, [3] , with the following expressions 

= 2 exp/3 (D + 2KQ)smh (2/3 JM) 

1 + 2exp[3(D + 2KQ) cosh (2/3 JM) 1 ' 



and 



or 



2 exp (3 (D + 2KQ) cosh (2/3 JM) 
Q ~ 1 + 2 exp /3(D + 2KQ) cosh (2{3JM) ^ ' 



M = Q tanh (2(3 J M) (6) 

where (5 = tj^. 

It was observed, from the analysis of these expressions, that, as a function of the S ratio, there are 
possible dipolar, quadrupolar or successive transitions. Assuming that K favors condensed state and 
J determines the ferromagnetic or antiferromagnetic ordering, we consider that this type of analysis 
is important for the identification of the microstructural properties which are responsible for the 
magnetic behavior of magnetic fluids. 

For a more detailed analysis of the properties of such a system, we represent the Bragg- Williams 
free energy / site in the form 

$ = -DQ - KQ 2 - JM 2 - kT {In A/" - (1 - Q) ln[A/"(l - Q)]- 

- 1(q + M )h\m(q + m)]-\(q-m)h\m(q-m)\^ (7) 

where N represents the number of sites of the lattice, and, for its equilibrium value, determined by 
the condition 



2 



is obtained the condition 



D - 2KQ - 2 J , /v 



l-^Q(i-Q) l 



or 




lfT Q(l-Q) Q + M 



2 j§M(l-Q) Q-M 



(9) 



fi 



ln2(l - Q) - In VO 2 - M 2 - - 2Q(1 Q) In - 
v ^' V ^ 2M(1-Q) 



'Q + M 
Q-M 



where the following notations are used 

2J _ 

kf ~ 7 



A /i 



(10) 



It was interesting for us to represent the curves D = D(Q), for various values of T, in different 
physical situations, determined by 7-values. 

At the first stage, from figure |l| we observe that, as a consequence of the mathematical (kinematical) 
relationship between Q and M, the only relevant 7-values are in the range 1.1 < 7 < 15. It is because, 
for 7 < 1.1, the Q-order parameter values less than unity are not founded, and, on the other hand, 
for 7 > 15, the curves are 7- insensitive, being almost of the same form. 
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Figure 1: Q = Q{M) dependence for various 7 values 

Our analysis reveals the 7— insensivity of the inflexion points for each D = D(Q) which preserves 
\x in the domain (—1, 0). Taking [i in the same domain, the analytical form of the curves is almost 
insensible to variations of \x. 

Comparing the figures that show de dependency D = D(Q), we can conclude that, for a given 
temperature, the form of the curves (extremes and inflexion points) hence, the existence of different 
types of phase transitions, is strongly influenced by the coupling constant D, while the J-coupling 
constant determines the degree of order at which the transition occurs (the increase of Q, and in the 
same time that of M). 
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It is obviously that the restriction for H = excludes the discussion about the segregation phe- 
nomena in ferrofluids. The problem will be analyzed in a future paper. 

Considering the structural complexity of magnetic fluids we proposed, |l(J, regarding phase tran- 
sitions, the use of an Ising-3/2 type model. For that we assumed the existence of "free" ferroparticles, 
i.e. ferroparticles which, by means of various causes, are not completely covered with surfactant. The 
concentration of this "free" ferroparticles strongly influences the agglomeration process. 

In defining the model we assume the same cell construction, as used by Kalikmanov, Q, but 
with the needed specifications for a lattice model. The above considerations imply the existence of 
two types of cells: cells within we have ferroparticles covered with a surfactant layer and cells which 
contain "free" ferroparticles. Two spin states will be allocated for each type of cell: one state for 
the ferromagnetic particle with or without it's surfactant layer and the second state for the remaining 
solvent in the cell. The main reason of this construction lies in the different dimensions of ferroparticles 
and solvent molecules (the first are much greater than the latter). 

In the general case, for the Ising-S model, the general form of the Hamiltonian is, [11|, 



H 



'j 



(Si 



ls(s + i] 



s 4 



\s(s + i) 



(12) 



The Ising-3/2 Hamiltonian can now be obtained now considering S = 4 in (12) 



n 



-jj^SiSj-KY: 



v 



(S l 



(13) 



In the case S = | it is also possible to introduce an "octupolar" ordering, with an order parameter 
given by (Si) 3 , but this is not the case here, since the Hamiltonian contains no terms of the form 
(Si) 3 (Sj) 3 . If we introduce the molecular fields H and D, associated with the two order parameters 
M and Q, the one-particle Hamiltonian in the mean field approximation will be 



Ho 



(S> 



(14) 



As it was shown, the Ising-3/2 Hamiltonian can be rewritten so that, the equations for M and Q 



in the mean field approximation, became [11 



M 



Q 



3e (3(2KQ+D) sin ^(S/3JM + l(3H) + e P(~2KQ+D) sinh ^J M + lpjf) 

2e p(2KQ+D) cos h(3/3JM + p#) + 2e Pi-^Q+D) cosh(/3JM + \[5H) 

2e 2f3(2KQ+D) cosh ( 3/3 J M + 3p H }_ 2e (3(-2KQ+D) cosh ^J M + 1^) 

2e l3(2KQ+D) cos h(3/3JM + Mh) + 2e^~ 2K Q +D ) cosh(/?JM + \(3H) 



(15) 
(16) 



As a result, we get a greater number of situations in which phase transitions occur. 

Same considerations as in the spin-1 case, in the absence of an external magnetic field, led us to 
the conclusion that the Q = Q(M) dependence is strongly determined by the value of D, and, in 
some cases, we get two values for the order parameter Q, both corresponding to the same value of 
M. As an example, for a given set of coupling constants in the domain of interest (7 = 1, \i = 0), 
the Q = Q(M) dependence has the form given in figure ^, from which we observe that, in addition to 
the M-discontinuity domain 0.74 -7- 0.8, at Q = 0, it also exists a nonanalitical behavior of M in the 
domain of configurational ordering with Q / 0. 
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Figure 2: Q = Q{M) dependence 
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